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Motivated by recent findings of unconventional superconductors exhibiting multiple broken symmetries, we 
consider a general Hamiltonian describing coexistence of itinerant ferromagnetism, spin-orbit coupling and 
mixed spin-singlet/triplet superconducting pairing in the context of mean-field theory. The Hamiltonian is di- 
agonalized and exact eigenvalues are obtained, thus allowing us to write down the coupled gap equations for 
the different order parameters. Our results may then be applied to any model describing coexistence of any 
combination of these three phenomena. As a specific application of our results, we consider tunneling between 
a normal metal and a noncentrosymmetric superconductor with mixed singlet and triplet gaps. The conductance 
spectrum reveals information about these gaps in addition to how the influence of spin-orbit coupling is mani- 
fested. Explicitly, we find well-pronounced peaks and bumps in the spectrum at voltages corresponding to the 
sum and the difference of the magnitude of the singlet and triplet components. Our results may thus be helpful 
in determining the relative sizes of the singlet and triplet gaps in noncentrosymmetric superconductors. We also 
consider the coexistence of itinerant ferromagnetism and triplet superconductivity as a model for recently dis- 
covered ferromagnetic superconductors. The coupled gap equations are solved self-consistently, and we study 
the conditions necessary to obtain the coexistent regime of ferromagnetism and superconductivity. Analytical 
expressions are presented for the order parameters, and we provide an analysis of the free energy to identify 
the preferred system state. It is found that the uniform coexistence of ferromagnetism and superconductivity 
is energetically favored compared to both the purely ferromagnetic state and the unitary superconducting state 
with zero magnetization. Moreover, we make specific predictions concerning the heat capacity for a ferromag- 
netic superconductor. In particular, we report a nonuniversal relative jump in the specific heat, depending on 
the magnetization of the system, at the uppermost superconducting phase transition. We propose that this may 
be exploited to obtain information about both the superconducting pairing symmetry realized in ferromagnetic 
superconductors in addition to the magnitude of the exchange splitting between majority and minority spin 
bands. 

PACS numbers: 74.20.-z, 74.25.-q, 74.45.+C, 74.50.+r, 74.20.Rp 



I. INTRODUCTION 

Recent findings of superconductors that simultaneously 
exhibit multiple spontaneously broken symmetries, such as 
ferromagnetic order or lack of an inversion center''^'"* and 
even combinations of such broken symmetries'*, have led to 
much theoretical and experimental researchSsM. The sym- 
metry of the superconducting gap in these and other uncon- 
ventional superconductors is presently a matter of intense 
investigation*^ '"'"''^''^^. Multiple spontaneously broken sym- 
metries are not only of interest in terms of studying prop- 
erties of specific condensed matter systems, but also due to 
the fact that it may provide clues for what could be expected 
in other systems in vastly different areas of physics. Topics 
such as mass-differences of elementary particles and emer- 
gent phenomena in biology is caused by spontaneously bro- 
ken symmetries''*, and in many cases, the phenomena may 
even be described by the same type of equations. In this pa- 
per, we will address the issue of competition and coexistence 
between three phenomena giving rise to broken symmetries 
which are highly relevant in condensed-matter physics: feiTO- 
magnetism, superconductivity, and spin-orbit coupling. 

The discovery of superconducting materials that lack a 
centre of inversion^"*'^"^ '^, such as CePtaSi, Ulr, Li2Pd3B, 
Li2Pt3B, and Cd2Re207, has lately triggered extensive the- 
oretical work on these compunds. Properties of a supercon- 



ductor without an inversion center were investigated early by 
Edelstein'*, while in Ref. [l^ it was shown that a 2D su- 
perconducting system with a significant spin-orbit coupling 
induced by the lack of inversion symmetry would display a 
mixed singlet-triplet superconducting state. This means that 
the superconducting order parameter would possess the ex- 
otic feature of having no definite parity. Later studies '^-'^-^° 
also investigated specific noncentrosymmetric superconduc- 
tors with a model Hamiltonian consisting of a superposition 
of spin-orbit and superconducting terms. In an attempt to de- 
termine the coiTect pairing symmetry of the superconducting 
state in such unconventional superconductors, it was found 
that the favored triplet pairing stata^ for the heavy-fermion 
material CePtaSi is dk oc {ky, —k^, 0). Very recently, how- 
ever, an experimental study^^ of thermal transport properties 
in the present compound concluded that the correct gap func- 
tion (dk-vector) may exhibit nodal lines in contrast to the 
point nodes displayed by the dk-vector suggested by Ref.l2Tl 
It is therefore of considerable interest to investigate several 
specific models for noncentrosymmetric superconductors in 
order to reveal characteristic features in physical observables 
that might be helpful in classifying the symmetry of the su- 
perconducting order parameter. 

In Ref. im the authors studied tunneling between a normal 
metal and a noncentrosymmetric superconductor considering 
the particular form of dk suggested by Ref. ,21, in the limit 



2 



of weak spin-orbit coupling and in the absence of spin-singlet 
pairing. Anderson^"* showed that the only stable triplet pair- 
ing states in the presence of a spin-orbit coupling would have 
to satisfy dk || gk, where gk = — g-k is the vector func- 
tion describing this interaction, such that in CePtsSi one also 
has gk = ^{ky, —kx, 0). Moreover, it was demonstrated by 
Samokhin^'' that the spin-orbit coupling in this particular ma- 
terial is significant, i.e. ^ ksTc, which indicates admixtur- 
ing of singlet and triplet Cooper pairs. In the present paper, 
we solve the full Bogoliubov-de Gennes (BdG)-equations for 
a system with spin-orbit coupling including both spin-singlet 
and spin-triplet superconducting gaps, studying a gap vector 
dk, point {ky, —kx,0) as suggested by Ref. i21i. We then ap- 
ply this gap vector to what we believe is the simplest model 
that captures the essential features that could be expected to 
appear in the conductance spectrum of a 2D normal/CePtaSi 
junction. Our work then significantly extends the considera- 
tions made in Ref. ^2? primarily in that we present analytical 
and numerical results that allow for both triplet and singlet 
gap components. Also note that a similar Hamiltonian was 
very recently studied in Ref. 25, where it was shown that the 
presence of a weak external magnetic field would significantly 
change the nodal topology of CePtsSi. With regard to non- 
centrosymmetricity, we underline that breaking the symmetry 
of spatial inversion does not in general give rise to a signif- 
icant spin-orbit coupling. Also, it is well-known that spin- 
orbit coupling may be induced in a centrosymmetric crystal 
by means of an external symmetry-breaking electrical field. 
In the latter case, however, the broken symmetry is strictly 
speaking not spontaneous as it certainly is for e.g. a crystal 
lattice undergoing a structural phase transition which breaks 
spatial inversion"^. 

Another interesting scenario in the context of sponta- 
neously broken symmetries is the study of superconductors 
that exhibit coexistence of ferromagnetic and superconduct- 
ing order, i.e. systems where two continuous internal sym- 
metries SU{2) and U{1) are simultaneously broken. Due to 
the preferred orientation of the spins in a ferromagnetic sys- 
tem, the SU{2) rotational symmetry is spontaneously broken. 
In a superconducting system, the ground state spontaneously 
breaks the U{1) symmetry. Note that by the terminology bro- 
ken symmetry, we are referring to the fact that the wavefunc- 
tion describing the state of the system acquires a complex 
phase which characterizes the ground state. In the ferromag- 
netic and superconducting systems we will consider in this 
paper, superconductivity appears at a lower temperature than 
the temperature at which onset of ferromagnetism is found. 
This may be simply due to the fact that the energy scales for 
the two phenomena are quite different, with the exchange en- 
ergy naturally being the largest. It may, however, also be due 
to the fact that superconductivity is dependent on ferromag- 
netism for its very existence. Such a suggestion has recently 
been put forth^^. 

In the context of FMSCs, it is crucial to address the question 
of whether the superconductivity and ferromagnetism order 
parameters coexist uniformly or if they are phase-separated. 
One plausible scenaric^^ is that a spontaneously formed vor- 
tex lattice due to the internal magnetization M is realized. 



but studies of a uniform superconducting phase in spin-triplet 
FMSCs^ has also been conducted. As argued by Mineev in 
Ref. [3^ an important factor with respect to whether a vor- 
tice lattice appears or not should be the magnitude of the in- 
ternal magnetization M. Specifically, Ref. 31 suggested that 
vortices may arise if 47rM > Hd, where Hd is the lower 
critical field. In the case of URhGe, a weakly ferromagnetic 
state coexisting with superconductivity seems to be realized, 
and the domain structure in the absence of an external field 
is thus vortex-free. Unfortunately, current experimental data 
concerning URhGe are not as of yet strong enough to unam- 
biguously settle this question. On the other hand, evidence for 
uniform coexistence of ferromagnetism and superconductiv- 
ity has been indicated^^ in UGe2. 

Although this is an unsettled issue, it seems natural to as- 
sume that in ferromagnetic superconductors (FMSCs), the 
electrons involved in the SU{2) symmetry breaking also 
participate in the U{1) symmetry breaking. As a conse- 
quence, uniform coexistence of spin-singlet superconductiv- 
ity and ferromagnetism can be discarded since s-wave Cooper 
pairs carry a total spin of zero, although spatially mod- 
ulated order parameters could allow for magnetic s-wave 
superconductors'^-^'*. However, spin-triplet Cooper pairs are 
in principle perfectly compatible with ferromagnetic order 
since they can carry a net magnetic moment. There is strong 
reason to believe that the correct pairing symmetries in the dis- 
covered FMSCs constitute non-unitary states'^''"'. Spin-triplet 
superconductors have a multicomponent order parameter dk, 
which for a given spin basis reads 



Ik = 



Akii-AktT -i(Aki + AkTT) 



A 



(1) 



Note that dk transforms like a vector under spin rotations. 
The superconducting order parameter is characterized as non- 
unitary if i(dk X dj.) 7^ 0, which effectively means that time- 
reversal symmetry is broken in the spin part of the Cooper 
pairs, since the average spin of Cooper pairs is given as 
(Sk) = i(dk X dj^). Notice that time-reversal symmetry 
may be broken in the orbital part (angular momentum) of the 
Cooper pair wavefunction even if the state is unitary. In the 
general case where all SC gaps are included, it is generally 
argued that Ak|x would be suppressed in the presence of a 
Zeeman-splitting between the f, j conduction bands. Distin- 
guishing between unitary and non-unitary states in FMSCs is 
clearly one of the primary objectives in terms of identifying 
the correct order parameter. Studies of quantum transport in 
junctions involving FMSCs has explicitly shown that the con- 
ductance spectrum should be helpful in revealing the correct 
pairing symmetry'^'^^. Hence, an itinerant electron model of 
ferromagnetism augmented by a suitable pairing kernel should 
be a reasonable starting point for describing such systems. 

Although we have mentioned two specific examples of sys- 
tems exhibiting multiple broken symmetries, our aim with this 
paper is to construct a solid starting point for consideration 
of a condensed-matter system exhibiting any combination of 
the broken symmetries resulting from superconductivity, fer- 
romagnetism, and/or spin-orbit coupling. By applying the ap- 
propriate limits to our theory, one may then obtain special 
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cases such as FMSCs or noncentrosymmetric superconduc- 
tors with significant spin-orbit coupling. 

This paper is organized as follows. In Sec. HIl we establish 
the Hamiltonian accounting for general coexistence of ferro- 
magnetism, spin-orbit coupling, and superconductivity. The 
diagonalization procedure and coupled gap equations are de- 
scribed in Sec. |III| Then, we apply our findings to a model 
of normal/noncentrosymmetric superconductor junction, cal- 
culating the tunneling conductance spectrum Sec. IIVI in addi- 
tion to a discussion of these results. As a second application, 
we consider a FMSC in Sec. [V] solving the coupled gap equa- 
tions self-consistently and calculating the free energy and heat 
capacity of such a system. Our main conclusions are summa- 
rized in Sec. |Vll We will use boldface notation for vectors, .?. 
for operators, and .T. for 2x2 matrices. 



II. MODEL FOR COEXISTENCE OF 
FERROMAGNETISM, SPIN-ORBIT COUPLING, AND 
SUPERCONDUCTIVITY 

For our model, we will write down a Hamiltonian describ- 
ing the kinetic energy, exchange energy, spin-orbit coupling, 
and attractive electron-electron interaction, respectively. The 
total Hamiltonian can then be written as 



respect to spatial inversion symmetry read 



where the respective individual terms read 

^kin = ^ Ckcj^-^Cktj, 

i7FM = -JAf^7(k)Sk-S_k, 

k 

^^SOC = ^ 4a (gk ■ Cr)al3C\^p, 
ka/3 

ffsC = (14Va,3 + ^kk'a/3)cLc-k/3C-k'/3Ck' 



(2) 



kk'ct/3 



(3) 



Above, Ck = £k — M where ek is the dispersion relation for 
the free fermions and /i is the chemical potential^^, J > 
is a ferromagnetic coupling parameter, 7(k) is a geometrical 
structure factor for the lattice, gk is a vector function account- 
ing for the antisymmetric spin-orbit coupling, while 14;k'Q/3 
is an attractive pair potential. The factor of 1/2 in iJsc is 
included to obtain more convenient expressions later on, and 
simply corresponds to a redefinition of V^^,^p 
In Eqs. (O, the spin operators are given by 



lyS.T 

2 "^kk'Q/J- 



Sk = 



N ^ 



i,3C(k+k')/3- 



(4) 



Moreover, we have explicitly split the attractive pairing po- 
tential into a singlet and triplet part according to Vkk'a/g = 
^kk'a/3 + ^kk'a/3- ^he Symmetry properties of the antisym- 
metric spin-orbit coupling and superconductivity terms with 



gk = -g-k, 



kk'Q/3 



= V: 



±k±'k'a/3j 



^kk'a/3 — ± ±' V'^k±'k'a/3- 



(5) 



In order to find eigenvalues and gap equations for our sys- 
tem, we introduce the mean-field approximation for the two- 
particle Hamiltonians (ferromagnetic and superconducting 
terms) such that the operators Sk and C]{j.^cLk^ may be writ- 
ten as a mean-field value pluss small fluctuations. We define 

(4acLk/3) = 4a,3' ^d writC 



(Sk) 



<5(Sk) 

^t 



(6) 



Inserting Eqs. (|6]l into Eqs. (|3]l and discarding all terms of 
order 0{6'^), one obtains in the standard fashion 



INAP 



ka/3 



ka,3 

+ [^Lp + ^LfMJ-up-bLp}- (7) 

In Eqs. (|7]l, M = {M^,My,Mz) = (Sj) = (Sk=o) denotes 
the mean value of the spin operators in real space, interpreted 
as the magnetization of the system. We have introduced the 
vector describing the magnetic exchange energy Vm — M 
and the order parameters (OPs) 

V = (Vm). - i(VA/)y = - iMy), V, = IM,, (8) 

for ferromagnetism, while the OP for superconductivity is de- 
scribed by 

1 



A S.T _ J_ T^S.T 7 

k' 



(9) 



The quantity / appearing in Eq. (HJ is a measure of the 
strength of the magnetic exchange coupling. Although we 
have derived the ferromagnetic part of our Hamiltonian from 
a lattice-model [where / = 2 JiV7(0)], this generic Hamilto- 
nian describes a general mean-field model of a system with 
magnetic exchange energy. The Pauli principle places the fol- 
lowing restrictions upon the superconductivity OPs: 



Singlet pairing: A^^ = -A^^„, A^^ = Aika/3- 
Triplet pairing: A^^ = A^^„, A£„^ = -Alka/3- 



(10) 



In total, we have thus obtained a Hamiltonian H de- 
scribing coexistence of ferromagnetism, spin-orbit coupling, 
and superconductivity in the mean-field approximation by 
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adding all of the above terms. For more compact nota- where we have introduced the quantities 
tion, one may introduce a basis for fermion operators (/)k = 
[ckT, Cki, cl^p cL^J'^ and write 

H = -ffkin + HpM + Hsoc + Hsc 

= ^O + ^^0l^k0k, (11) 
k 



Ho 



k 

V(AT 



CkT + 5k, ^ 
-Vt + .gk, 

(A^tt)^ 



kTi 



^ y^(Akaj3 + AkQ/3)&kQ/37 



ka/3 



Cki 



.9k.- 

5k,2 



(AjTl)^-(A^Ti) 



^kTT 
^kTi 



A^ 
^kTT 



A^ 

^kTi 



(A^k 



kii^ 



^CkT + .9k, z 

V + .9k,- 



+ .9k,+ 

-Cki - gk,z 1 



(12) 



Above, we have defined Cka = Ck — crVz in addition to 
gk,± = (gk)a; ± i(gk)y The matrix Ak will be central in 
this work, and we note that it may be further compactified by 
introducing the dk-vector formalism^. By means of the defi- 
nitions c?k,o = Akij^ and 

dk = i[A^^T - Ai^Tp-i(AT^^ + AT^^),2AT^J, (13) 

that transforms like a vector under spin rotations, one may 
write 

A ^ /^Cki - • (Vm - gk) i(rfo.k + dk • 0-)(7y \ 

UiKk + dk-(T)<7j,]t -Cki + (Va/ + gk) • <tV ' 

(14) 

where i denotes the identity matrix and T designates the ma- 
trix transpose. The rest of this paper wiU now be devoted 
to obtaining the excitation energies for H by diagonalizing 
Ak, writing down the coupled gap equations, and considering 
some important special cases. 



III. EXCITATION ENERGIES AND GAP EQUATIONS 

The characteristic polynomial for a general matrix with 
eigenvalues iJk may be written as^ 

+ - /| + |Ak + /|] - 1 - det^kji^^ 

+ - /| - l^k + 1\\E^, + detAk = 0, (15) 

where / denotes the 4x4 identity matrix. Since in our case 
is Hermitian, Tr{Ak} = 0, and the polynomial reduces to a 
depressed quartic equation. For ease of notation, we introduce 



the quantity 

^k ^^[l^k-/|±|^k + /|], (16) 

such that Eq. (fTsT i is rewritten as 

K + {^l - 1 - detAk)^;^ + 77k -Sk + detAk = 0. (17) 
The solutions of i?k can be written as^ 

2E^^^ = aak + /3{-[3(r/+ - 1 - det^) 

+ 2yk + «^]r/'. (18) 
ak 

Here, we have defined the auxiliary quantities 

ak = V^l^yk - 1 - detAk) + 2yk, 

5(r/+ - 1 - detAk) , 

2/k = — -„ - bk, 

6 

h - pl/^ R - "Sk , , / Qk I -Pk 

(77+ - 1 - detAk)detAk (7?k)' 

3 8- 

_ (7?+ - 1 - detAk)' 
108 

(^1^ - 1 - detAk)^ 
Pk = — - detAk. (19) 

In Eq. (fTsl l. {a, /?} take the values +1 and —1 such that there 
exists a total of four solutions for i?k- Also note that any of 
the roots in the expressions for 6k and i?k will do the job. A 
special case of the above solutions, which occurs quite fre- 
quently in various contexts, considerably simplifies the ob- 
tained eigenvalues: ^[lAk — | — |^k + -^|] = 0. In this case. 



5 



the quartic equation reduces to an effective quadratic equation 
with the solutions 



- 2(r/+ - 1 - detAk) 



2/?V ('/k - 1 - detAk)2 - 4detAk 



1/2 



(20) 



This is the situation considered in most problems dealing with 
superconductors. Having calculated the energy eigenvalues, 
Eq. ( fTTT i may now be diagonalized by writing 



H = Ho 

= Hq 



k 

jE('^k^k)(P^^kPk)(P^0k) 



k^k^k. 



(21) 



where Z?k = diag(i?k.i; -Ek,2, -Ek.a, i^k^) is a diagonal ma- 
trix containing the eigenvalues of A\^. Here, we have defined 
[see Eq. ([H] 



1 1 



(22) 



thus absorbing the factor ^ in front of into the eigenval- 
ues. Above, Pk are the orthonormal diagonalizing matrices 
which by the hermiticity of ^k are ensured to be unitary. We 
write our new basis of fermion operators as 



t t 

•^k = [7kT'7ki'7-kT'7-ki] 



(23) 



These operators satisfy the fermion anticommutation rela- 
tions, as can be verified by direct insertion. From Eq. ( |2TI) . 
we may now write 

H ^Ho + EbkT7kT-E^k,i + 7k;7ki£^k.2 

k 

+ 7-kT7lkT-^k,3 + 7-ki7-ki-E^k,4] 

= [Ho + E(^k,3 + £^k,4)] 
k 

+ E[^kT'^kT(i^ka-S-k,3) 

k 

+ 7,^x7ki(£^k,2-^-k,4)] 
= Ha + Et'^kT^kT-Bk,! + 7ki7ki^k,2], (24) 

k 

where we have defined Hq — Hq + X]k(^k,3 + -E'k.4) and 

Ew,i = (£;ka-i^-k,3),i^k,2 = (£^k.2-S-k,4). OurHamil- 
tonian now has the form of a free-fermion theory. It is then 
readily seen that the free energy of the system is given by 

F^Ho-^ EP"(1 + c^^^'^ O + ln(l + c-''^'^-^)]. (25) 



From F, the gap equations for the ferromagnetic and super- 
conducting OPs V,Vz, and A^'^^ may be obtained by de- 
manding the value of these which corresponds to a minimum 
in F. The possible extrema of F are given by the conditions 



dF 



dA 

OF 



S,T 
Waf3 

= 0, 



0, 



dF 



9(Ak«/3r 
dF OF 



0, 



= 0. 



(26) 



By first defining the quantity 



dE 



'k,l 



dEi 



'k,3 



+ nF(£^k,2) 



dx 

dEk,2 
dx 



dx 

dx ^ 



(27) 



where nf{E) — 1/(1 + e'^-^) is the Fermi distribution, the 
conditions in Eqs. ( |26] | may be evaluated by inserting Eq. 
dZSl l. The extrema of F are thus determined by the following 
equations: 



— bka/S 

2J7(0) 

N 

-Y + 

N 



0, 



4- ^(V.) 
NV 



= 0, 

T{V) = 0, 



0. 



(28) 
(29) 

(30) 

(31) 
(32) 



The challenge then lies in obtaining the derivatives of the ener- 
gies i?k,i with respect to the different order parameters. In the 
general case described by Eq. (fTSl) . this is a formidable task. 
Nevertheless, the above above provides a general framework 
which may serve as a starting point for any model considering 
the coexistence of ferromagnetism, spin-orbit coupling, and 
superconductivity. We will apply our findings onto a specific 
case which currently is a topic attracting much attention: non- 
centrosymmetric superconductors with significant spin-orbit 
coupling. 



IV. PROBING THE PAIRING SYMMETRY OF 
NONCENTROSYMMETRIC SUPERCONDUCTORS 

As an application of our model, we consider tunneling 
between a normal metal and a noncentrosymmetric super- 
conductor treated in the spin-generahzed Blonder- Tinkham- 
Klapwijk (BTK) formalisi"^— 



A. Model and formulation 

The Hamiltonian in the superconducting state using stan- 
dard mean-field theory with a spin-orbit coupling may be writ- 
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ten as 



1 



k 



(33) 



using a spin basis <^k = [ckji Ck]., cL,^-., c^j^, and with 



Mk = 



3k,+ 
^kTT 



TT 



.9k- Ak 
fk -Ak 
-A|, -ek 



9K- 



Akii 
5k, + 



(34) 



In Eq. ( |34] |. all quantities have been defined in the previous 
section. It is usually argued that interband pairing in a noncen- 
trosymmetric superconductors can be neglected due to a spin- 
split Fermi surface in the presence of spin-orbit coupling. This 
is motivated by realizing that the splitting could be as large 
as^^ 50-200 meV for the noncentrosymmetric superconduc- 
tor CePtaSi, thus far greater than the superconducting critical 
temperature /csTc — 0.06 meV in that compound. Accord- 
ingly, one might be tempted to also exclude the spin-singlet 
gap in the presence of a strong spin-orbit coupling motivated 
on physical grounds by the suppression of interband-pairing 
due to the spin-split Fermi surfaces. However, it is necessary 
to investigate the presence, although possibly small in mag- 
nitude, of a spin-singlet component of the gap to examine 
whether the conductance spectrum changes significantly in 
any respect compared to the scenario with exclusively triplet 
pairing. Another motivation for including the singlet gap is 
that the authors of Ref . 2 1 demonstrated that for small spin- 
orbit coupling, dk II gk yields the highest Tc for CePtaSi. 
This would thus correspond to a scenario where the triplet 
gap Ak|x is suppressed due to the above condition, although 
intraband-pairing is not strictly forbidden as a result of weak 
spin-orbit coupling, thus allowing for singlet pairing. 

Consider now a gap vector exhibiting point nodes. Since 
dk in general is given by Eq. ([T]), the vector characterizing 



spin-orbit coupling gk = A(fcj,, —k^^ 0) suggested by Ref.l2ll 
results in 

A, At 

^^^^ " ^21^*-^^ ^ ^^^^ " 2|k[''^^ ~ 

Diagonalization of the Hamiltonian in Eq. ( |33] ) yields eigen- 
values and eigenvectors which are necessary to calculate the 
normal- and Andreev -reflection coefficients in a N/CePtsSi 
junction. Assuming the simplest form of a s-wave supercon- 
ducting gap that obeys the symmetry requirements dictated 
by the Pauli-principle, namely an isotropic gap Ak = A,, we 
find that the eigenvalues of Afk read 



£^k„, - av/(£ + /3|gk|)2 + |A, + /3A,/2|' 



(36) 



This is in complete agreement with the result of Ref. |25l We 
are here assuming that all gaps have the same phase associated 
with the broken C/(l) gauge symmetry. In Eq. ( l36b . a = 
+ (— ) refers to electronlike (holelike) excitations, while j3 = 
+ (— ) denotes the spin-orbit helicity index. The wavevectors 
may then be written as 



qI = qJi — \/ kp + m^A^ — mX, 



i = qi = J kp + m^A^ + mX, 



(37) 



when making the approximation that the magnitude of the su- 
perconducting gaps is small compared to the Fermi energy /i 
and considering the low-energy transport regime. Here, kp is 
the Fermi wave-vector 

We now calculate the normal- and Andreev-reflection coef- 
ficients for an incident electron with spin a, which in turn will 
allow us to derive the tunneling conductance of the junction. 
To do so, we first set up the Bogoliubov-de Gennes (BdG)- 
equations for the system which read (see Appendix A for a 
derivation): 



' 2m 



X{py 



fi + VoSix) x{py + ipx)e{x) 
-Ate(x) 



AkTT0(a;) 
-A,e(a;) 



Akiie(x) 



A[^^e(a;) 



+ fl-VoSix) X{Py - 1P^J;)Q{X) 



Ale{x) 



Al,e{x) 



X{py + iPx)0{x) 



*(a;,y) = E^{x,y), 



(38) 



M - VoSix)/ 



where Px(y) = ^'^^^iy) make use of the boundary condi- 
tions 

;) V(0) = ^"(0) (Continuity of wavefunction), 

ii) 2mVo^{0) = 4*(0) - 4V^'(0) - mA77^'(0) 
(Continuity of flux). (39) 

Note that we have applied the usual step-function approxi- 
mation for the order parameters instead of solving for their 



spatial dependence self-consistently near the interface, i.e. 
Ako-o-(r) = Ako-o-6(a;) and A(r) = A8(x) (we comment 
further on this later). For convenience, we have defined the 
4x4 matrix 



1 0\ 
-10 
1 
0-10/ 



(40) 



7 



The presence of spin-orbit coupling leads to off-diagonal 
components in the velocity operator, such that it would be er- 
roneous to merely match the derivatives of the wavefunction 
in this case"*^. The coupled gap equations that arise by de- 
manding a minimum in the free energy are obtained by con- 
sidering Eqs. dZTl i and ( |28] |. For the sake of obtaining ana- 
lytical results, we continue our discussion of the conductance 
spectra of noncentrosymmetric superconductors by inserting 
values of the superconductivity gaps a priori instead of using 
the self-consistent solutions. This approach does not, then, 
account for the entire physical picture, but has proven to yield 
satisfactory results for many aspects of quasiparticle tunneling 
in the case of e.g. spin-singlet d-wave superconductors^'''"*'*. 

For the simplest model that illustrate the new physics, we 
have thus chosen a two-dimensional N/CePtaSi junction with 
a barrier modelled by V{y) — Vq6{x) and superconductiv- 
ity gaps Akcrcr(r) = Akcrcr0(a;) [5{x) and 6(x) represent 
the Delta- and Heaviside-function, respectively]. Consider 
Fig. [T] for an overview. Choosing a plane- wave solution 
'^{x, y) = 5'(x)e''"'''^, for a =t the wavefunction on the nor- 
mal ['(/'(a;)] side of the junction reads 



il){x) = 



/giAiF COS 9x 

V 



7.ip — i^F COS 6x 
o^^'F COS Ox 
r^ikp COS Ox 



1 



(41) 



On the superconducting side [^'(2:)], the BdG-equation may 
be written, for our particular choice of gk and gaps in Eq. 
\, as 



/ 




l5k|e'^ 


-(A,/2)e'* 


As \ 




l<7k|e-'* 


£k 


-A, 


(At/2)e-'* 




(A,/2)e-'^ 






l5k|e-'^ 


V 


A, 


(A./2)e'* 


l.9k|e'^ 


-Ek / 



= with tan 0(6*) = l/tan6'. 

(42) 



We are here concerned with positive excitations E > 0, as- 
suming an incident electron above Fermi level. In this case, 
there are four possible solutions for wavevectors k with a 
given energy E > 0. Consequently, one may verify that the 
correct wavefunction for x > 0, which is a linear combination 
of these allowed states, reads 



V2 



/ «(A+) \ 
u(A+)e-'"*(''') 
-w(A+)e-"^(^=) 

V '^(A+) J 
( v{A+) \ 

w(A+)e-"^(<) 
-u(A+)e-'"^(<) 

V u{A+) J 



lijj cos ej X I 

V2 



T T f 
iq; COS e,[ X _| h_ 

^/2 



/ u(A_) \ 
^ -m(A_)c-'^(«') ' 
v{A^)e-"t'^o^^ 
\ «(A-) / 
/ <A.) \ 
-i;(A_)e-"^(''h) 
u(A_)e-"^(«h) 
V u{A.) J 



gigj" cos 9^x 



^ig,! cos 9lx 



(43) 



We have defined A± = As± ( At/2), and the spreading angles 
in Eq. (l43T l are given as sin6'g = kp sinO / q" , 6'^ = tt — 
9'^. This follows from the fact that translational symmetry is 
conserved along the y-axis. The coherence factors entering 
the wavefunctions in Eq. ( |43] | are given as 

.(A) ^ ^/iT^^, .(A) . ./iT^^ 
^ ^ V 2 2E ' ^ ' \ 2 2E 

(44) 

We also define the dimensionless parameters Z = 2mVo/kp 
and (3 — 2mX/kp as a measure of the intrinsic barrier strength 
and magnitude of the spin-orbit coupling, respectively. 

Note that we are using the same effective masses in the nor- 
mal part of the system as in the superconducting part. The 
mass of the quasiparticles in heavy-fermion materials are, as 



the name itself implies, ordinarily much larger than in normal 
metals. It was recently shown by Yokoyama et al.*^ that in 
a two-dimensional electron gas (2DEG)/superconductor junc- 
tion where spin-orbit coupling was substantial in the 2DEG, 
the effect of including a larger effective mass in the 2DEG was 
equivalent to that caused by an increase of Z. Note that in 
the presence of a time-reversal breaking magnetic field, it was 
shown in Refs. 46 47 that the effect of Fermi-vector mismatch 
could not be reproduced simply by varying the barrier param- 
eter Z. Since there is no time-reversal breaking field present 
in this case, however, we here restrict ourselves to consider- 
ing equal effective masses in the two systems. With the above 
equations, one is able to find explicit expressions for {r^ , r^}. 
The procedure illustrated here is identical for incoming elec- 
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Normal metcil 



Reflected electron (T, J.) 



Incoming electron a 



Retroreflected hole (T, i) 




Noncentrosymmetric superconductor 

Transmitted e-likc qnasiparticlc {a) 
Traiismitted e-like quasiparticle (—a) 

Ttaiismitted e-like quasiparticle (-cr) 



Transmitted e-Kkc quasiparticle (a) 



FIG. 1: (Color online) Schematic illustration of the scattering pro- 
cesses taking place at the interface of the 2D planar N/CePtsSi junc- 
tion. The arrows indicate the direction of group velocity (which is not 
equal to the momentum vector for the holes). Note that the presence 
of spin-orbit coupling causes electron-like and hole-like excitations 
on the superconducting side to be spread into different angles. 



trons with a —I when using 

/ 

e 

V 



ikp cos Ox 



-ikp cos 6x 



cjl^F COS 6x 

oifcp COS 6x 



(45) 



instead of Eq. ( |4TI ). This establishes the framework which 
serves as the basis for calculating the conductance spectrum. 



B. Conductance spectra for noncentrosymmetric 
superconductors 

We now proceed to calculate the tunneling conductance for 
our setup. Generalizing the theory of Blonder, Tinkham, and 
Klapjwik"^^, one obtains a conductance G{E, 9) (scaled on the 
conductance in a N-N junction) for an incoming electron with 
angle 9 to the junction normal with spin a, where 

G{E, 9) = l + Y,{\rl{E, 9)\^ KiE, 9)W (46) 

a 

and Rn-n = ^^^^^^[A^cos^ 6'/(4cos^ 9 
larly averaged conductance reads 



Z'^)]d9. The angu- 



G{E) = {Rn-n) 



G{E,9)P{9)d9, (47) 



where P{9) is the probability distribution function [P(0) = 1] 
for incoming electrons at an angle 9. This is in many cases set 
to P{9) = cos 9, but other forms modelling e.g. effective 
tunneling cones may also be applied. In obtaining the total 
conductance, one has to find G{E) for both a and a —[ 
and then add these contributions. The original derivation of 
this sgecific formula for the tunneling conductance given in 
Ref.WareUes on the relation 



to hold. This is known to be valid for subgap energies, but 
for energies above the gap the relationship does not hold in 
general, a fact which implies that the conductance formula 
derived in Ref. 48 is only valid for applied voltages below 
the gap, strictly speaking. However, since the probability for 
Andreev reflection rapidly diminishes for energies above the 
gap (especially for Z ^ 0), the conductance formula may still 
be applied for larger voltages as a reasonable approximation, 
even for the high-transparency case of low values for Z. 
The explicit analytical expressions for the normal- and 



Andreev-reflection probabilities. 



and 



respectively. 



are too large and unwieldy to be of any instructive use. We 
shall therefore be content with plotting these expressions to re- 
veal the physics embedded within them. In most scanning tun- 
neling microscopy (STM) experiments, a high transparancy 
interface is often realized, corresponding to low Z. Also, 
since the band-splitting 2Afci? at Fermi level may be of order— 
100 meV for CePtaSi, a simple analysis relating this to our di- 
mensionless parameter (3 yields that (3 ~ 0.05. We therefore 
plot in Fig. |2]the angularly averaged (and normalized) con- 
ductance spectrum for several values of barrier strength and 
singlet/triplet gap ratios, fixing the spin-orbit coupling param- 
eter ?A j3 — 0.05. From Fig. |2] we see that one may infer the 
relative size of the singlet and triplet components of the gap by 
the characteristic behaviour of G{E) at voltages correspond- 
ing to As ± At/2. This is in agreement with what one could 
expect by studying the form of the eigenvalues in Eq. ( l36b . 
since it is this precise combination of the gaps that appear in 
the expression. 

In a recent study '^'^ by Iniotakis et ai, a nor- 
mal/noncentrosymmetric superconductor junction was 
studied for low-transparency interfaces, where it was found 
that zero-bias anomalies would take place for certain STM 
measurement orientations if a specific form of the mixed 
singlet-triplet order parameter was realized. This may be 
attributed to the formation of zero energy bound states^, 
which is possible when the gap contain nodes. In the present 
study, we are using an isotropic spin-singlet gap and and also 
isotropic p-wave gaps (|Akcrcr| = constant)^^, such one does 
not expect the appearance of a ZBCP, in contrast to Ref. '4^ 
Moreover, we note that the spin-orbit coupling in the system 
gives rise to effectively spin-active boundary conditions [see 
Eq. dllll]^. 

It is also instructive to consider the Andreev-reflection 
probabilities explicitly to resolve the spin-structure of the 
quasiparticle current, as shown in Fig. [3] for an incoming 
electron with spin a It is seen that the spin-| coefficient 
becomes larger with increasing voltage, such that the spin- 
polarization of the current will vary with the bias voltage. The 
proper definition of a spin-current in systems exhibiting spin- 
orbit coupling has, however, been shown^^ to be more subtle 
than applying the usual relations for charge- and spin-currents 



J charge ^ 



(49) 



(48) 



where is the particle-current of fermions with spin a. 
Therefore, it is fair to claim that it is not obvious how one 
might detect such a change in polarization of the quasiparticle 
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FIG. 2: (Color online) Tunneling conductance for N/CePtaSi junc- 
tion with l3 = 0.05. We study barrier strengths corresponding to a) 
Z = 0.1, h) Z — 1, c) Z — 10. In all cases, we plot the ratios 
At/ As = {3, 10} to see how the spectra are affected. It is seen that 
the conductance spectra reveal information about the relative size of 
the singlet and triplet component of the gaps by characteristic fea- 
tures located at bias voltages £ = |As ± At/2j. 



current with a change in bias voltage. On the other hand, the 
charge-current remains unaffected by these considerations and 
our results thus indicate that the conductance spectrum of the 
charge-current in a N/CePtaSi junction may provide valuable 
information about the relative size of the singlet and triplet 
components of the superconductivity gap. 

We now comment on effects that have not been taken into 
account in our analysis of this problem. First, the issue of 
how boundary effects affect the order parameters is addressed. 
Studies^^i^i^ have shown that interfaces/surfaces may have a 
pair-breaking effect on unconventional superconductivity or- 
der parameters. This is relevant in tunneling junction exper- 
iments as in the present case. The suppression of the order 
parameter is caused by a formation of so-called midgap sur- 
face states (also known as zero-energy states)^^ which occurs 
for certain orientations of the k-dependent superconducting 
gaps that satisfy a resonance condition. Note that this is not 
the case for conventional s-wave superconductors since the 
gap is isotropic in that case. This pair-breaking surface ef- 
fect was studied specifically for p-wave order parameters in 
Refs.'53!'54', and it was found that the component of the order 
parameter that experiences a sign change under the transfor- 
mation k± —k±, where fcj^ is the component of momen- 
tum perpendicular to the tunneling junction, was suppressed 
in the vicinity of the junction. By vicinity of the junction, 
we here mean a distance comparable to the coherence length, 
typically of order 1-10 nm. Thus, depending on the explicit 
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FIG. 3: (Color online) Andreev-reflection coefficients for spin-| and 
spin-| fermions in the case of incoming a =| electrons. It is seen 
that the degree of spin-polarization of the generated quasiparticle 
current will vary with the bias voltage. The inset contains a plot of 
the sum of reflection coefficients (both normal and Andreev), show- 
ing that no transmittance of quasiparticles occurs for voltages below 
|As- A,/2|. 



form of the superconducting gaps in the noncentrosymmetric 
superconductor, these could be suppressed close to the inter- 
face. Moreover, we are dealing with an easily observable ef- 
fect, since distinguishing between the peaks occuring for vari- 
ous values of Ra requires a resolution of order 0(10^^ A|_o)> 
which typically corresponds to 0.1 — 1 meV. These structures 
should readily be resolved with present-day STM technology. 
However, it should be pointed out that a challenge with re- 
spect to tunneling junctions is dealing with non-idealities at 
the interface which may affect the conductance spectrum. 

In order to fully consider the possible pair-breaking effect 
of the interface in an enhanced model, one would obviously 
need to solve the scattering problem self-consistently in order 
to obtain more precise results for the conductance, especially 
in terms of the quantitative aspect. To obtain analytical re- 
sults, however, we have inserted the gaps a priori, since we 
believe that our model captures essential qualitative features 
in a N/CePtaSi junction that could be probed for This belief 
is motivated by studies^** for (1^2 _y2 superconductors which 
show that the conductance shape around zero bias remains 
essentially unchanged even if the spatial dependence of the 
order parameters are taken into account. The spectra around 
the gap edges may be modified in the sense that since the gap 
in general will be somewhat reduced close the interface, the 
appearance of characteristic features in the conductance could 
occur at lower bias voltages than the bulk value of the gaps. 
However, it seems reasonable to hope that our simple model 
may be of use in predicting qualitative features of the conduc- 
tance spectrum when considering junctions involving noncen- 
trosymmetric superconductors such as CePtaSi. 
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V. PROBING THE PAIRING SYMMETRY OF 
FERROMAGNETIC SUPERCONDUCTORS 

As a second application of our model, we consider a model 
of a ferromagnetic superconductor described by uniformly co- 
existing itinerant ferromagnetism and equal-spin pairing non- 
unitary spin-triplet superconductivity. 



A. Model and formulation 

We write down a mean-field theory Hamiltonian with 
equal-spin pairing Cooper pairs and a finite magnetization 
along the easy-axis similar to the model studied in Refs. ,59i. 
leolleii namely 



k 2 2 ^ ^k<T<TOk<Ta 

k k<T 



^kcr 



A 



kcrcr 



^ko-o- ^kcr 



Ckcr 



(50) 



Applying the diagonalization procedure described in Sec. [Ill 
we arrive at 



k(T 

^0 = ^ - ^k. - AL.&k..) 



(51) 



k(T 



where {7kcri7kcr} new fermion operators and the eigen- 
values read 



k(T 



+ |Ak. 



(52) 



Recall that it is implicit in our notation that is measured 
from Fermi level. The free energy is obtained by using the 
procedure explained in Sec. [Ill and one obtains 



(53) 



k<T 



such that the gap equations for the magnetic and supercon- 
ducting order parameters become^- 



M = - 



N ^ 



kCT 



tanh(/3.Ek./2), 



kk'crcr 



2-Bk' 



-tanh(/3£k'./2). (54) 



For concreteness, we now consider a specific form of the gaps, 
similar to those studied in Refs. l59ll6ll Assuming that the gap 
is fixed on the Fermi surface in the weak-coupling limit, we 
write 



Akcrcr — Ar 



A^^o 



(55) 



where ki? is the normalized Fermi wave-vector, such that the 
gap only depends on the direction of the latter We have intro- 
duced the spherical harmonics 



Y, 



i=i\ 



-ay/S/Sne"''^ sine 



(56) 



such that the gaps in Eq. dSSl l experience a change in sign 
under inversion of momentum, i.e. 6 ^ 6 + tt. We shall 
consider the case sine/) 1 which renders the magnitude 
of the gaps to be constant, similar to the A2-phase in liq- 
uid ■^He. The motivation for this is that it seems plausible 
that uniform coexistence of ferromagnetic and superconduct- 
ing order may only be realized in thin-film structures where 
the Meissner (diamagnetic) response of the superconductor is 
suppressed for in-plane magnetic fields. This enables us to 
set sin0 = 1, since the electrons are restricted from moving 
in the z-direction in a thin-film structure. In a bulk structure, 
as considered in Ref. 61, we expect that a spontaneous vortex 
lattice should be the favored thermodynamical state^^. The 
pairing potential may then in general be written as 

m 

v,A0,0') ^-J^Jf^Y^ewie')]*, (57) 

^ — ' O/OTT 

m ' 

which for the chosen gaps reduces to 



(58) 



Conversion to integral gap equations is accomplished by 
means of the identity 

^E/(^k.)= / deN'^ie), (59) 

k •' 

where N"'{e) is the spin-resolved density of states. In three 
spatial dimensions, this may be calculated from the dispersion 
relation by using the formula 



V 



(27r)3 7ek„=const IVkEko 



(60) 



With the dispersion relation ^kcr = £k ^ crIM — Ep (having 
set the chemical potential equal to the Fermi energy, /i = Ep), 
one obtains 



N^ie) = mVV2He + aIM + Ep_)^ ^^^^ 



27r2 

In their integral form, the gap equations read 



de 



1 = 



tanh[/3^,(e)/2], 
(62) 



I / de—y-tanh[PE.{e)/2]. 

B. Zero temperature case 



Consider now T — Q, where we are able to obtain analytical 
expressions for the superconductivity order parameters in the 
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problem. Since the superconductivity gap equation reduces to the self-consistency equation for the magnetization becomes 



1 = ^ 
2 



de 



(63) 



one readily finds 



(64) 



where we have defined M = lAI/Ep, i-e. the exchange en- 
ergy scaled on the Fermi energy. Moreover, the weak coupling 
constant c = gN{{))/2 will be set to 0.2 throughout the rest of 
this paper, unless specifically stated otherwise. Moreover, we 
set (Do — i^ol Ep — 0.01 as the typical spectral width of the 
bosons responsible for the attractive pairing potential. From 
Eq. (l64l i. we see that the effect of increasing the magnetiza- 
tion is an increase in the gap for majority spin. The impor- 
tant influence of the magnetization is that it modifies the den- 
sity of states, which affects the superconductivity gaps. For 
M = 1, i.e. an exchange splitting equal to the Fermi energy, 
the minority spin gap is completely suppressed, as shown in 
Fig. |4] Thus, the presence of magnetization reduces the avail- 
able phase space for the minority spin Cooper pairs, suppress- 
ing the gap and the critical temperature compared to the pure 
Bardeen-Cooper-Schrieffer (BCS) case. 



x10 



5- 



Magnetization enhances the majority spin gap 



Zero DOS for minority spin at Fermi level 



' - - ^ A|,o 



0.2 



0.4 0.6 

M 



0.8 



FIG. 4: (Color online) Superconducting gaps (full drawn line: ma- 
jority spin, dashed line: minority spin) as a function of the magneti- 
zation with ujo = O.OI-Ef- When the exchange splitting equals the 
Fermi energy, the DOS of minority spin fermions is zero at Fermi 
level, resulting in a complete suppression of A^.o- 

After the appropriate algebraic manipulations of Eq. ( l62b . 



f{M) = M 



I 



1 - 29 - 



A2 (Af) 



1 - 26 



(x + 2A/)2 + A2p(A/) 




{x + 2M)-^J{x + 2MY + AJ o(Af) 



where we have defined the parameter / — IN{0), in simi- 
larity to Ref. 59, and introduced A^^o{M) = Aafl/Ep- We 
have thus managed to decouple the gap equations completely, 
such that one only has to solve Eq. i65[ to find the magnetiza- 
tion, and then plug that value into Eq. (l64b . Note that strictly 
speaking, one should divide the integral in Eq. ( l65T l into three 
parts: { — 1 — M,—ujo},{ljo,oo},{—ujo,ujo} where the su- 
perconductivity gaps are only non-zero in the latter interval. 
However, the error associated with doing the integration nu- 
merically over the entire regime with a finite value for the gaps 
is completely neglible. From Eq. ( |65] l, we see that the trivial 
solution AI = is always possible. Interestingly, we find that 
a non-trivial solution implying coexistence of ferromagnetism 
and superconductivity is only possible when / > 1 (in agree- 
ment with Ref. 59). To illustrate this fact, consider Fig. |5]for a 
plot of f{M) in Eq. ( |65] ) as a function of M for several values 
of /. In fact, it is seen that more than one solution is possi- 
ble for any / > 1: the trivial solution M = Q corresponding 
to a unitary superconducting state, and a non-trivial solution 
M — Afo, representing a non-unitary superconducting state. 
Recall that in terms of the dk-vector formalism, these classi- 
fications are defined as 

Unitary: dk x = 0, Non-unitary: dk x dk 7^ 0. (66) 

We will later show that the free energy is minimal in the non- 
unitary state, which implies that the coexistence of ferromag- 
netism and superconductivity may indeed be realized in our 
model. 

The order parameters depend on the parameters (T,I,c). 
To illustrate their dependence on / at T = 0, consider Fig. 
|6] It is clearly seen that the superconductivity gaps are equal 
for / < 1, corresponding to a unitary spin-triplet pairing 
state. For / > 1, a spontaneous magnetization arises and the 
majority/minority spin gap increases/decreases. This corre- 
sponds to the coexistent phase of ferromagnetism and super- 
conductivity. An important point concerning Eq. i65[ is the 
inclusion of the step-function factors, which are superfluous 
as long as we are considering the coexistent regime of fer- 
romagnetism and superconductivity, since their argument is 
always negative. However, if one for instance were to set one 
or both of the superconductivity gaps to zero, the correct gap 
equation for the magnetization would not be reproduced with- 
out them. This is due to the loss of generality in taking the 
limit tanh{(3Ecr) 1 when /3 — > 00 in deriving Eq. ( |65] ). 
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X 10" 




an unconventional superconductors. As we have alluded to, 
it seems reasonable to assume that a superconducting phase 
analogous to the Al- or A2-phase of "^He may be realized in 
FMSCs. We now investigate how the magnetization at T = 
depends on the ferromagnetic exchange energy constant / in 
these possible phases, and compare them to the purely ferro- 
magnetic case. Our results are shown in Fig. [T] where we 
have self-consistently solved for M as a function of / in three 
cases: 1) the purely ferromagnetic phase, 2) the Al-phase 
where only spin-f fermions are paired, and 3) the j42-phase 
where all spin-bands participate in the superconducting pair- 
ing. It is seen that the magnetization is practically identical 
in all phases regardless of the value of /. Analytically, this 
may be understood since the difference A/ [see Eq. (l65T l1 
between the gap equation for the magnetization in the purely 
ferromagnetic case and the coexistent state reads 



FIG. 5: (Color online) Plot of f{M) in Eq. ([gill as a function of 
M. Only the trivial solution exists for 7 < 1, while three possible 
solutions are present at / > 1. We have plotted Eq. l l65b for values 
of [0.9,1.1] in steps of 0.01. 



since Ecr > is replaced with e when superconductivity is 
lost, which can be both larger and smaller than zero when 
Act.o 0. The present form of Eq. $65[ is generally valid 
for the purely magnetic and the coexistent Al- and A2-phases 
of the ferromagnetic superconductor 



X 10" 



deiV' 



A2 



0. 

(67) 



Note that in our results, an enhancement of the magnetiza- 
tion below the superconductivity critical temperature is ab- 
sent, contrary to the results of Ref. |6lj who predicted that the 
magnetization should be enhanced in the coexistent phases 
compared to the purely feiTomagnetic phase. For the weak- 
coupling approach applied here, it seems reasonable that the 
presence of superconductivity should not alter the magnetiza- 
tion much, while superconductivity itself is drastically mod- 
ified depending on the strength of the exchange energy. The 
result of Ref. 61 may be a consequence of the fact that they do 
not set sin = 1 [Eq. (|56]ll, and consequently have additional 
nodes compared to the gaps we are using. 




0.99 



FIG. 6: (Color online) Self-consistently solved order parameters at 
T = as a function of I. It is seen that the coexistent regime of 
ferromagnetism and superconductivity is located at J > 1, where a 
spontaneous magnetization arises. 

In order to correctly characterize the pairing symmetry of 
FMSCs, it is of interest to find clear-cut experimental sig- 
natures that distinguish between the possible phases of such 



1.035 
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FIG. 7: (Color online) Plot of the exchange energy M = IM/Ef 
at r = as a function of J = IN{0) for two possible phases of 
a FMSC (the Al- and yl2-phase) compared to the purely ferromag- 
netic case. It is seen that M is virtually unaltered by the presence 
of superconductivity, at least in the weak-coupling approach we have 
adopted here (see also Ref.l6ll). 
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C. Finite temperature case 

The critical temperature for the superconductivity order pa- 
rameter is found by solving the equation 

1 = 1/ de ^tanh £/2T,,, , (68) 

which yields the BCS-like solution 



T,,, = l.UiVoe-'/'^^+^^^^K (69) 

Since the transition temperature for paramagnetism - ferro- 
magnetism is in general much larger than the superconduct- 
ing phase transition, one may to good approximation set 
M{Tc^cr) = M{0). It is then evident that the critical tem- 
perature depends on the magnetization in the same manner as 
the gap itself, and the cutoff-dependence in Eq. ( l64b may be 
removed in favor of the critical temperature by substituting 
Eq. ( |69] l. In order to solve the coupled gap equations self- 
consistently at arbitrary temperature, we considered Eq. ( l62T i 
with the result given in Fig. [8] It is seen that the minority- 
spin gap is clearly suppressed compared to the majority-spin 
gap in the presence of a net magnetization. Also, the graph 
clearly shows that the BCS-temperature dependence consti- 
tutes an excellent approximation for the decrease of the OPs 
with temperature. In what follows, we shall therefore use self- 
consistently obtained solutions at T = for the OPs and make 
use of the BCS temperature-dependence unless specifically 
stated otherwise. In general, the critical temperature for the 
ferromagnetic order parameter, Tc,m exceeds the supercon- 
ducting phase transition temperatures Tc,a by several orders 
of magnitude. However, for / very close to one, we are able to 
make these transition temperatures comparable in magnitude. 
In the experimentally discovered FMSCs UGe2 and URhGe, 
one finds that Tc,m is 50-100 times higher than the tempera- 
ture at which superconductivity arises. 

To illustrate how the magnetic order parameter depends on 
/, consider Fig. |9]for a plot of the temperature dependence for 
several values of I. The inset shows how the critical tempera- 
ture depends on this parameter. 



D. Comparison of free energies 

Although a non-trivial solution of M exists, care must be 
exercised before concluding that this is the preferred energet- 
ical configuration of the system. Specifically, it may in the- 
ory be possible that the systems prefers the M — solution 
regardless of the value of /, corresponding to a unitary super- 
conducting state with A|,o = '^i.o- It is therefore necessary 
to compare the free energies of the M = and M ^ cases 
at values of / where the latter is a possible solution, and also 
study their temperature dependence. In the general case, the 
analytical expression for the free energy in the coexistent non- 



X 10"" 



/= 1.01 

lO-^M 




10^ 10"' 10"" 10^ 10"^ 10"' 



T 

FIG. 8: (Color online) Self-consistently solved order parameters as 
a function of temperature for I = 1.01. Note that the tempera- 
ture axis is logarithmic, such that the transition between the para- 
magnetic and ferromagnetic phase is much higher than the super- 
conducting phase transitions. However, we are able to tune / such 
that Tc,M and T^.a become comparable. We have also plotted the 
gaps with self-consistently solved values at T = and then apply- 
ing a BCS-temperature dependence (solid black lines), which yield 
excellent consistency with the solution that does not assume a BCS- 
temperature dependence. 




FIG. 9: (Color online) Temperature-dependence of the magnetic or- 
der parameter for several values of /. The insets shows how the 
critical temperature depends on /. 
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unitary superconducting phase reads 



F/N 



„ A9 „ J-Ep-alM 



^' + ^1, 1 



(70) 



Note that the gap should be set to zero in the above equation 
everywhere except in the interval {— wq, ^o}- We obtain a di- 
mensionless measure of the free energy by multiplying with 
I/Ej,, and denote i^NU = FI/{NEp). Note that the free en- 
ergies of the unitary state, pure ferromagnetic state, and para- 
magnetic state are obtained as follows: 

Fu = lim Fnu, 

Af^O 



F 



FM — lini -Fnu- 



In Fig. [To] we plot the difference between the unitary and non- 
unitary solution at zero temperature, AF = F\j ~ Fnu, which 
clearly shows how the system favors the non-unitary solution 
with spontaneous magnetization as / increases. As a result, 
we suggest that the coexistent phase of ferromagnetism and 
superconductivity should be realized at sufficiently low tem- 
peratures whenever a magnetic exchange energy is present. 
For consistency, we also verified that Fnu < ^fm at T = 
since the system otherwise would prefer to leave supercon- 
ductivity out of the picture and stay purely ferromagnetic. 
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FIG. 10: (Color online) Comparison between the free energy for the 
non-unitary and unitary superconducting state at zero temperature. It 
is seen that these values are equal for 7=1 (and for / < 1), while 
the non-unitary state is energetically favored for increasing ferromag- 
netic exchange energy. Thus, the coexistent phase should be realized 
at sufficiently low temperatures in the presence of a ferromagnetic 
exchange energy. 

We now turn to the temperature-dependence of the free 
energy at the fixed value of / = 1.01 (the order parame- 
ters were self-consistently solved for this value and plotted 



in Fig. [8]l. The results are shown in Figs. [TT]to[T3] Note 
that we now use a different scaling of the free energy, namely 
Fnv = F/[NN{0)Tl^]. The well-known result that the free 
energy of a purely superconducting state joins the free energy 
of the paramagnetic state continuously as the temperature in- 
creases is reproduced in Fig. [TT] In Fig. [121 we see that 
the coexistent phase of ferromagnetism and superconductiv- 
ity is energetically favored compared to the purely ferromag- 
netic case, which is consistent with the experimental fact that a 
transition to superconductivity occurs below the Curie temper- 
ature for certain materials'-^. Finally, in Fig. [13] we have plot- 
ted the energy difference between the unitary and non-unitary 
free energy in addition to the difference between the paramag- 
netic and ferromagnetic phases. It is seen that the non-unitary 
state is energetically preferred over the unitary state, a state- 
ment which strictly speaking has only been shown to hold for 
our current choice of / (/ = 1.01), but it seems reasonable 
to assume that it holds under quite general circumstances due 
(71) 

to the presence of an exchange energy. At T — Tcj, when 



all superconductivity is lost, the two curves join each other 
smoothly since Fnu ^fm and Fu FpM when T > Tc,]. 
Our results then suggest the very real possibility of a coexis- 
tent phase of spin-triplet superconducting pairing and itiner- 
ant ferromagnetism being realized in the experimentally dis- 
covered ferromagnetic superconductors, since we have shown 
that the coexistent phase is energetically favored over both the 
purely magnetic and the non-magnetic superconducting state. 



0.14 




FIG. 1 1 : (Color online) Free energy difference between the param- 
agnetic state (Fvm) and the unitary superconducting state (Fu). In 
consistency with established results (see e.g. Ref. 62), the free en- 
ergies merge continuously as the temperature gets closer to Tc.u. In 
the inset, we have chosen the zero-temperature value of the param- 
agnetic free energy as zero, serving as a reference point. We have 
solved all order parameters self-consistently for / = 1.01. 



E. Specific iieat 

We next consider some experimental signatures that could 
be expected in the different possible phases of a FMSC. Con- 
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FIG. 12: (Color online) Free energy difference between the fer- 
romagnetic state (Ffu) and the non-unitary superconducting state 
(Fnu) which displays coexistence of ferromagnetism and supercon- 
ductivity. It is seen that the non-unitary phase is favored compared 
to the purely ferromagnetic state. In the inset, we have chosen the 
zero-temperature value of the ferromagnetic free energy as zero, 
serving as a reference point. We have solved all order parameters 
self-consistently for / — 1.01. The curves of Figs. 12 and 13 may 
be made congruent by a simple scaling of the axes. This is a conse- 
quence of the weak-coupling limit, where superconductivity sets in at 
a temperature much smaller than the ferromagnet-paramagnet transi- 
tion temperature, such the that magnetic order parameter across the 
superconducting transition essentially is a temperature-independent 
constant. 
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FIG. 13: (Color online) Free energy difference between the unitary 
and non-unitary state (Fu — Fnu) as well as the paramagnetic and 
ferromagnetic state (Fpm — Jfm). At T = 7c, the curves merge 
smoothly into each other since all superconductivity is lost. Each 
step along the ordinate corresponds to an increment of unit 0. 1 . We 
have solved all order parameters self-consistently for / = 1.01. 



sequently, we have calculated the electronic contribution to 



the specific heat of the system by making use of Cy — T 



dS 
dT 



with 

S = - Y.if(E^-Mf{Eu.)] + [l-f{Eua)]H^^f{Eua)]} 

(72) 

as the entropy, leading to 



^ EL 



dT 



cosh^(/3Sk./2) 



(73) 



Note that the above equation reduces to the correct normal- 
state heat capacity in the limit {Ao-.c-A^} 0, with the 
usual linear T-dependence. The term 



dT 



ensures that 



the well-known mean-field BCS discontinuity (strictly speak- 
ing valid only for a type-I superconductor^"*, but clearly 
invalid at the transition temperature of a strong type II 
superconductor''"'-^^-^^) at the superconducting critical temper- 
ature is present in the heat capacity, while the presence of fer- 



dM 
dT ■ 



How- 



romagnetism induces a new term proportional to 
ever, due to our previous argument that Tc^m ^ Tc.cr, this 
term may be neglected since the magnetization remains virtu- 
ally unaltered in the temperature regime around Tc ^. Going 
to the integral representation of the equation for the heat ca- 
pacity, one thus obtains 



deN%e)[El{e)--^/^„,,T] 
cosh-2[/3^,(£)/2]. (74) 



Strictly speaking, one should again divide the above inte- 
gral into the regions {—Ep, — wq}, {woi oo}, and {— ajg, t^o] 
where the superconductivity gap should be set to zero in all 
regions except the latter. However, since the integrand is 
strongly peaked around e = (Fermi level), there is little 
error made in using the form Eq. (f74l l. In order to obtain the 
derivatives of the gap functions with respect to temperature, 
an analytical approach is permissable since the gaps have the 
BCS-temperature dependence (see Fig. |8]l 



A„,o(r) = A,,o(0)tanh 1.74. /T.^./T - 1 



(75) 



where the superconductivity critical temperature for spin-cr 
fermions is given by Eq. ( |69] ). To illustrate how the super- 
conductivity pairing symmetry leaves important fingerprints 
in the heat capacity, we solved Eq. ( |74l i self-consistently for 
two values of / coiTesponding to strong (Af ~ 0.5) and weak 
{M ~ 0.1) exchange splitting. At / = 1.01, the disconti- 
nuity is clearly pronounced for T — Tc.], but it is hardly 
discernable at T = T^^i- However, for / — 1.0005 where 
the superconductivity transition temperatures for majority and 
minority spins become comparable, a clear double-peak sig- 
nature is revealed in the heat capacity. We thus propose that 
this particular feature should serve as unambigous evidence 
of a superconducting pairing coiTesponding to the A2-phase 
of liquid '^He in ferromagnetic superconductors. 

An classic feature of the BCS-theory of superconductivity 
was the prediction that the jump in the heat capacity at Tc 
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FIG. 14: (Color online) Specific heat capacity as a function of tem- 
perature for two values of /, corresponding to a strong exchange 
splitting (A/ ~ 0.5) and a weak exchange splitting (A/ ~ 0.1). 
A double-peak signature is clearly visible when the transition tem- 
peratures for the majority and minority spin bands are comparable. 



normalized on the normal-state value was a universal number, 
namely 



/ ACy 
V Cv 



1.43. 



(76) 



In the presence of a net magnetization, one would expect that 
the universality of this ratio would break down and depend 
on the strength of the exchange energy. This is due to the 
fact that the discontinuity in the specific heat at the supercon- 
ducting transition is dominated by the majority-spin carriers, 
while the total specific heat to a larger extent has contribu- 
tions from both minority-spin and majority spin carriers. To 
investigate this statement quantitatively, we consider the jump 
in Cv at r = r^T since no analytical approach is possible 
at T = Tc_j, as seen from Eq. ( |74l l. We find that the normal 
(ferromagnetic) state heat capacity reads 



/-fFM _ 
Oy — 



7V-(0), 



(77) 



where A^'^(O) is the spin-resolved DOS at Fermi level, while 
the difference between the heat capacity in the coexistent state 
and the ferromagnetic state at T = Tc_| reads 



ACv 



1.742 A2(o)ArT(o) 



2T, 



Since the zero-temperature value for the gap is A 
l.TGTcT' arrives at 



V Cv J 



T = Tr 



= 1.43— 

■T 1 



1 



1-M 
1+Af 



(78) 

.0(0) = 

(79) 



The above equation reduces to the BCS-limit for complete 
spin-polarization M = 1 (zero DOS for spin-| fermions at 



Fermi level). This is due to, as noted above, the larger extent 
to which majority-spin carriers dominate the jump in specific 
heat compared to the total specific heat. As anticipated, the 
jump in Cv depends on the exchange energy, as illustrated in 
Fig. [15] Of course, in the unitary state M = the jump also 
reduces to the BCS value although this is not seen from Eq. 
(I79I 1. The reason for this is that we have implicitly assumed 
that M 7^ in the derivation of Eq. (|29]l, taking TcT > T^. 
In the case where these transition temperatures are equal, the 
contribution from both is additive and equal [ 1 .43/2, to be spe- 
cific, as seen from Eq. ( |79] ll and gives the correct BCS result. 
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FIG. 15: (Color online) The discontinuity of the heat capacity at 
T = Tc, I as a function of exchange splitting [Eq. J79H. It is seen 
that the BCS value is recovered at A/ = 1. Note that it would also 
be recovered at Af — 0, although this is not shown explicitely in 
the figure. The reason for this is that we have assumed that Tc,i 7^ 
Tcj. We have also plotted the numerical results (o) for the jump with 
self-consistently solved OPs, i.e. without assuming BCS temperature 
dependence, for I = {1.001, 1.005, 1.01, 1.02, 1.05}, which yield 
good agreement with the analytical solution Eq. ( |79b . 

Our study of Cv then offers two interesting opportunities: 
i) the presence or absence of a double-peak signature in the 
heat capacity reveals information about the superconductivity 
pairing symmetry realized in the FMSC, and ii) the normal- 
ized value of the discontinuous jump at Tc,! contains infor- 
mation about the exchange splitting between the majority and 
minority spin carrier bands. 



VI. SUMMARY 

In summary, we have derived a general Hamiltonian 
describing coexistence of itinerant ferromagnetism, spin-orbit 
coupling and mixed spin-singlet/triplet superconducting 
pairing using mean-field theory. Exact eigenvalues and 
coupled gap equations for the different order parameters have 
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been obtained. Our results may serve as a starting point for 
any model describing coexistence of any combination of 
these three phenomena simply by applying the appropriate 
Umit. 

As a specific application of our results, we have studied 
quantum transport between a normal metal and a supercon- 
ductor lacking an inversion center with mixed singlet and 
triplet gaps. We find that there are pronounced peaks and 
bumps in the conductance spectrum at voltages corresponding 
to the sum and difference of the magnitude of the singlet and 
triplet gaps. Consequently, our results may be helpful in ob- 
taining information about the size of the relative contribution 
of different pairing symmetries. 

Moreover, we considered a system where itinerant 
ferromagnetism uniformly coexists with spin-triplet super- 
conductivity as a second application of our theory. We 
solved the coupled gap equations numerically, and presented 
analytical expressions for the order parameters and their 
dependences on quantities such as exchange energy and 
temperature. It was found that the coexistent regime of fer- 
romagnetism and superconductivity may indeed be realized 
since it is energetically favored compared to a unitary super- 
conducting state (M = 0) and a purely ferromagnetic state. 
In order to make contact with the experimental situation, we 
studied the heat capacity and found interesting signatures in 
the spectrum that may be used in order to obtain information 
about both the superconductivity pairing symmetry present in 
the system and the magnitude of the exchange energy. 
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APPENDIX A: BOGOLIUBOV-DE GENNES EQUATIONS 
FOR SYSTEMS EXHIBITING COEXISTENCE OF 
FERROMAGNETISM, SPIN-ORBIT COUPLING, AND 
SUPERCONDUCTIVITY 



1. Derivation 



We start out with a real-space Hamiltonian described by 
fermionic field operators tp'^''\r,t) with a general attractive 



pairing kernel V^/3(|r — r'|), namely 



a/3 



2m 



H + {V + aVs)S{x) 



+ [-VM + g(p)]-<^ i^l{r,t) 

J a/3 

+ J drdr'V^p{\r-r'\)i.i{r,t)^pl{r',t) 



X ^^(r',i)V'a(r,i). 



(Al) 



Here, Vq accounts for a non-magnetic scattering potential as- 
sociated with a barrier located at x = while Vg is the mag- 
netic scattering potential, i.e. the barrier is spin-active. More- 
over, Vm is the magnetic exchange energy vector, g(p) = 
— g(— p) is a term describing an antisymmetric spin-orbit cou- 
pUng energy (p = — iVr), while & is the vector of Pauli ma- 
trices. We now introduce the mean-field approximation 



(A2) 



where the last term describes the flucuations around the aver- 
age field, and also define the superconducting order parameter 



A„/3 {r,v') = V^p{\r-r'\) {^jp (r' , t) (r , t)) . 



(A3) 



Above, we have exphcitly made the superconductivity or- 
der parameter time-independent, which effectively amounts to 
saying that it does not depend on energy (the weak-coupUng 
limit). This provides us with 



H 



a/3 



t) 



2m 



-fx+{V + aV,)S{x) 



[- Vm + g(p)] • <t] Jl{r,t) 



In I /drdr'[At^(r,r')V/3(r',i)Va(r,t) 

af3 



-Vc0{\r-r'\){^i{r, t)^l (r' , t)) (V'/3 (r' , Va (r , i))] . 



(A4) 



The time-dependent field operators ip{r, t) = e 
obey the Heisenberg equations of motion 

idttlJa{r,t) = [ipa{r,t),H] 

+ f dr'A„/3(r,r')V^(r',t), 

= Y.f dr'^(r-r')[-^°(r',-p)]J^vMr',t) 
+ E|dr'At^(r,r')V/3(r',t). (A5) 
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For convenience, we have defined 



iJ" (r,p)= --JL^^+{y + aVs)5{x) 



2m 

[-Vm + g(p)] • a- 



(A6) 



The above equations may be comprised in compact matrix 
form 

iat*(r,t) = j dv'n{v,v')-^{r',t), 

*(r,t) = [V't(i-,i),V^i(r,i),V^|(r,t),^|(r,t)]T, 



H(r,r') = 



A(r,r') 



-A (r,r') hi/"(r',-p)]^^, 



with (5(r — r') = 6rr', and where we have defined 



lAxT(r,r') Au(r,r') 



(A7) 



(A8) 



Note that A|^(r, r') is in general a superposition of a triplet 
(T) and singlet (S) component that satisfy 

ATi(r,r') = AT^(r,r') + A^^(r,r'), 
A^,(r,r')-Aj^(r,r'), 



A?,(r,r') = -AS (r,r') 



(A9) 



Regarding ^ (r, t) as a c-number and assuming a stationary 
solution ^(r, i) = ^'(r)e~'^* with i? as the wavefunction 
energy, it suffices to solve the equation 

£;*(r) y"dr'7i(r,r')*(r'). (AlO) 

By considering a plane-wave solution of 5'(r) and dividing 
out the fast oscillations on an atomic-scale (see e.g. Ref. l63l) . 
one is left with most familiar form of the BdG-equations ap- 
pearing in the literature, namely 

fH°(r,p) A(k,r) \ , , 

At, r iror ',irPW=£^*r, (All) 
\A (k,r) [-H°{r,-p)Yj 

where the quasiparticle momentum k is the Fourier-transform 
of the relative-coordinate s = (r — r')/2, i.e. 

Hfm - / ds/(s)e-"'-. (A12) 



This is usually assumed to be fixed on the Fermi surface, such 
that only the directional dependence of k enters in Eq. dAl II ). 
k A;i?k. 



2. Boundary conditions 

We proceed to provide a general approach in order to obtain 
the correct boundary conditions at the interface for the wave- 
functions. Continuity of the wavefunction itself is assumed 
in this context. Consider our Eq. (l38l l which describes the 
Hamiltonian for the N/CePtaSi junction. The first row of the 
equation explicitly reads 



1 1 v6 ] 

d d 

+ - i—)Q{x)4'i {x, y) + AkTTe(a;)V'| {x, y) 



+ A,e(-x)^/.| (x, y) = Ei^^{x, y), 



(A13) 



If we now integrate the above equation over a an interval 
[e, — e] along the x-axis and apply the limit e 0+, one ob- 
tains 



lim 



{ - ^ [V-U^ ' y) - V-T (-e, y)] + i^oV-T (0- y) 

+ X f da;[e(x)^x(a;,y)]'} =0, (A14) 



where ' denotes derivation with respect to x. The last term 
yields ■iAV'x(e,0) (since 0(0) — i), such that the boundary 
condition for derivative of the V't {x, y) -component becomes 



^lim { (e, y) - {-e, y)] - m\^^ (e, 0) } 



= 2ml/oV'T(0-y)- 



(A15) 



It is seen that the presence of spin-orbit coupling and the delta- 
function barrier leads to a discontinuity of the derivative of 
the wave-function. A similar procedure may be applied to 
the other components of ^(x, y), and this method can also be 
extended to include different effective masses on each side of 
the junction modelled by a simple step-function 9(x). 
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